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1.0 INTRODUCTION 

The SAO work carried out in the period covered by this report consisted of 
a continuation of our effort to determine the size and shape of the hollow cathode 
cloud emitted from an orbiting system. In addition we have applied results 
obtained for the ionospheric wave impedance of a tethered system to the 
experiments under consideration. We are still considering the recent plasma 
chamber experimental results reported by Urrutia and Stenzel [1986(a) and (b)]. 
These UCLA investigators have been calling into question the operating principles 
and feasibility of electrodynamic tethered satellite systems, based on their results, 
which show currents to be limited to values lower than those predicted by probe 
theory. While they do not specifically consider a system with hollow cathode 
devices, their rather sweeping claims of a fundamental flaw in the electrodynamic 
tethered satellite system concept could reasonably be extended to include such 
systems. Perhaps these results will be a stimulus to performing some carefully 
designed, high magnetic field experiments on hollow cathode devices in a plasma 
chamber. There is a plasma chamber at JSC/NASA of about the same size as the 
one used in the UCLA experiments and already equipped with coils capable of 
giving an axial magnetic field of 0.4-40 G [J. McCoy, private conversation 1986]. 
Our final report will discuss the UCLA results in the context of designing plasma 
chamber experiments. 

One comment is worth making now. The supposed local circuit closure 
(across the magnetic field lines) observed in the UCLA experiments is a result that 
certainly cannot be extended to apply to the conditions of short tethers in space. 
The electrode separation in the UCLA plasma chamber was only 5 cm, or around 
30 electron gyroradii. For TSS-1 type conditions with a 200 m tether, the 
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separation would be — in terms of electron gyroradii — three orders of magnitude 
greater. 

The experimental results referred to above are not the only challenge raised 
recently against the ability of an electrodynamic tethered satellite system to draw 
currents above a few milliamps from the ionosphere. Barnett and Olbert [1986] of 
MIT analyzed the radiation from an electrodynamic tethered satellite system 
operating in a constant current mode. They found that wave impedances were 
far greater than previously estimated, reaching as high as 10,000-100,000 ohms. 
This was due to a hitherto neglected band of frequencies lying between the lower 
hybrid and electron cyclotron frequencies, which the motion of the system would 
supposedly excite. SAO has examined this problem and determined that the large 
wave impedances obtained by Barnett and Olbert do not apply to a real 
electrodynamic tethered satellite system but are a result of their having considered 
the unrealistic case of an orbiting wire. We briefly summarize our results in this 
report, since they are relevant to the experiments being considered. 


2.0 THE WAVE IMPEDANCE QUESTION 

In addition to the problem of local charge exchange across the sys- 
tem/plasma boundaries considered independently of the “circuit closure” there is 
the problem of how the ionospheric plasma’s global response to the moving 
disturbance represented by the electrodynamic tethered satellite system affects the 
system’s ability to function as a power generator or thruster. In the models of 
hollow cathode devices developed so far, no attempt has been made to take into 
account the fact that an electrodynamic tethered satellite system, even when the 
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current flowing through it is constant, causes a time-varying perturbation to the 
ionospheric plasma fields which generates electromagnetic plasma waves. A self- 
consistent theory that includes both the local plasma/system interactions and the 
large-scale perturbations caused by wave generation is a long way off, if not 
unattainable in practice. It is possible to estimate the effects of the plasma 
response on the functioning of the system, however, by calculating the wave 
impedance of the system under simplifying assumptions. 

Barnett and Olbert [1986] considered the problem of wave generation by an 
electrodynamic tethered satellite system and found wave impedances in the range 
of lOkfJ, mainly due to radiation in a frequency band lying between the lower 
hybrid and electron cyclotron frequencies. These unexpectedly large wave 
impedances, if correct, would restrict tether currents to well below an ampere, even 
with very long tethers. This limitation would apply to the experiments we are 
presently considering. 

SAO has been engaged in its own study of the tether wave generation 
problem, in mainly the context of deliberate, as opposed to incidental, generation 
of ULF/ELF waves. Our analysis can be found in detail in the final report for 
the NASA study NAG8-551. The SAO findings dispute the claim that an 
electrodynamic tethered satellite system would have a wave impedance of more 
than a few ohms. Since this issue is important for the short-tether experiments 
presently under consideration, we will summarize our results here. 

The main point that comes out of our analysis is that modeling well the 
dimensions of the charge-exchange interfaces between the system and the 
ionospheric plasma is of fundamental importance for correctly obtaining the 
electromagnetic wave fields excited in the plasma. The electrodynamic tethered 
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satellite system acts as a moving source (of opposite polarity at the two ends of 
the system) of external charge to the ionospheric plasma. It is the time-varying 
(because of the system’s motion) fields generated by this charge at the tether ends 
that excites the waves generated by the system. The plasma responds to the 
perturbation with electromagnetic waves which carry the net charge density away 
from the system. 

Thus wave generation is intimately connected with “circuit closure.” The 
field-line currents are electromagnetic plasma wave packets with nonzero diver- 
gence of the electric field. The frequencies associated with the tethered system’s 
perturbation of the ionosphere are determined by the time it takes for the charge- 
exchanging regions of the system to pass by the geomagnetic field lines. The 
larger the dimensions of these regions along the line-of -flight are, the lower is the 
maximum frequency associated with the disturbance. In practice these dimensions 
will be determined by the size of the terminating satellites or, more likely, by the 
size of the plasma contactor clouds. 

This conclusion, which might be reached directly by physical reasoning, also 
emerges from the equations used in our analysis. It is also in the equations of 
Barnett and Olbert, but they evidently didn’t draw the physical conclusion from it, 
since they modeled the tethered satellite system in an unrealistic way. They 
considered it to be a long, narrow cylinder, i.e., an orbiting wire. In effect they 
reduced the charge-exchange surface of the system to the cross-section of the wire 
at its two ends. It is the small dimensions of the charge-exchange region in their 
model that makes the higher frequency band so important in their analysis and 
causes them to obtain the very high wave impedances. For a system with charge- 
exchange region dimensions of a few meters, radiation in the band that is so 
important in the analysis of Barnett and Olbert becomes negligible. 
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This is not to deny that the Barnett and Olbert analysis is an important 
contribution. The exclusion of radiation (generated by a constant tether current) 
in the frequency band between the ion cyclotron and the lower hybrid frequencies is 
a significant conclusion (within the framework of cold plasma theory) that does not 
depend on the system dimensions used in the model. 

Wave impedances for a realistic model of the tethered system are around 
0.4fl for a tether 10 km long in the daytime central F-region. The dependence on 
tether length and terminating dimensions is weak. The wave impedance is linear 
in the Alfven speed, so increases by a factor of ten above the value quoted above 
might occur within a single orbit, based on the studies of electron density 
variation. However, the model is highly idealized. At this point its best use, 
perhaps, is to show that the Barnett and Olbert results, which were obtained with 
the same assumptions, are not applicable to a real tethered satellite system, and 
thus cannot be taken as proof that an electrodynamic tethered satellite system 
won’t work. 


3.0 SIMPLE PLASMA KINETIC MODEL FOR PLASMA CONTACTOR 
CLOUD* 

In our previous report for this investigation (Section 3.4) we considered 
the emission of the hollow cathode gas into the streaming ionospheric gas within 
the framework of a fluid model. This gave us a lower limit on the “standoff 
distance” to which the hollow cathode expands against the atmospheric stream. 
For the parameters considered in the earlier report this was only 20 cm. We note 
that this quantity varies inversely with the square root of the atmospheric mass 
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density. Increasing the altitude from the 220 km considered in the previous study 
up to 300 km (TSS-1 height) results in a ten-fold increase in the fluid model 
standoff distance to 2 meters. 

In the present analysis we estimate an upper limit to the standoff distance 
based on the opposite extreme of “weakly interacting” gas particles in the kinetic 
theory approach sketched below. 

Let the plasma contactor rest at the center of the coordinate system. We 
will work in a reference frame moving with the plasma contactor. 

Take the Shuttle flying in the + e x direction, and so the flow of ionospheric 
material past the plasma contactor cloud has velocity — V or b e x , where V or b is the 
orbital velocity of the Shuttle. 

Begin by writing down the ionospheric distribution function (in the absence 
of a plasma contactor cloud and the contactor cloud distribution function in the 
absence of interaction with the ionosphere). Thermal velocity spreads can be 
neglected for both distribution functions, at least for the initial treatment. 

The thermal spread of the ionospheric distribution function may be 
neglected because the Shuttle motion with respect to ionospheric material is highly 
supersonic. Consider a typical ionospheric species as represented by oxygen atoms 
at a temperature of 10 3 °K. 
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= 1.24 x 10 5 cm/ sec = 1.24 km / sec 

Vorb 8 km/sec _ 643 

Vth 1.24 km/ sec 

Although this is not an enormously large number, it is still probably possible 
to neglect ionospheric thermal velocities. 

The calculations on the adiabatic expansion of the plasma contactor cloud 
show the contactor cloud to be very cold, furthermore, its expansion velocity 

tfe<p w 2.5 x 10 4 cm / sec (see calculation in previous report) 


satisfies 


Vorb - 8xl0 5 cm/ sec _ 00 

^ — - — oZ 

Vex. 2.5 xlO 4 cm /sec 

Therefore, we may to an excellent approximation, neglect the thermal velocities 
of the plasma contactor cloud and also neglect any velocity dependence of collision 
cross-sections in treating particles traveling in the ± e* directions. 

Begin by adopting a model in which the interaction between the plasma 
contactor and the ionospheric flow is “weak.” In this model calculation wherever a 
scattering occurs between a contactor cloud particle and an ionospheric particle, the 
recoil velocities are large compared to any other velocities in the problem except 
the orbital velocity. The scattered particles are assumed to leave the system 
“instantaneously” without further scattering. This sort of “weak interaction 
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single-scattering” limit is clearly not realistic throughout the volume of the plasma 
contactor cloud and will yield an overestimate of the plasma contactor cloud size. 
This may be useful, however, as the fluid-dynamic calculation of the previous 
report bounds the size of the plasma cloud from below. It may be possible to 
obtain a reasonable estimate of the plasma contactor cloud bounded by these two 
limits. 


The distribution function of the ionospheric background, in the absence of 
interaction with the plasma contactor cloud may be written: 

F 6 (x,t/) = N 0 (x) 6(v x - Vorb ) 6 Ik) (1) 

It will be assumed that a steady state solution (d/dt = 0) can be found for the 
density distribution. Here N 0 (x) = background number density of particles 
(particles/cm 3 ), the spatial dependence of N 0 (x) (derivation from N 0 (x) = constant) 
is determined by scattering with the plasma contactor cloud. For purposes of this 
calculation it can be assumed that the background is composed of a single neutral 
atomic species, e.g. atomic oxygen. 

Define distribution functions for electrons, ions, and neutral atoms in the 
plasma contactor cloud as f e [x, 0), fi{x,v), and /„(x, t;) respectively. However, as 
ionization and recombination processes may be expected to be negligible outside the 
plasma contactor itself, we can refer to these generically as /(x, u). Again we are 
initially concerned only with stationary solutions. The (unperturbed) number 
density of the plasma contactor cloud particles is proportional to r -2 , owing to 
conservation of particle fluxes, i.e. 
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n(r) = no(r 0 ) (ro/r) 2 (2) 

where ito and r 0 axe reference values of number density and radius chosen to avoid 
the unphysical singularity at r = 0. The values of no and r<i are chosen on the 
basis of the mass flow rate and other properties of the plasma contactor, as 
outlined in the calculation of 10 October 1986. Then the plasma contactor cloud 
distribution function may be expressed in spherical coordinates as 

f{x,V) = n(r) S(Vr-Vexp) v eX p~ 2 (3) 

where v r is the radial velocity and v eX p the velocity of expansion of the plasma 
contactor cloud particles. Expressing the distribution function in rectangular 
coordinates 


f(z,v) — w(r) t)y -h tk ~ Vexp 'J/Vexp 

= n(r) + ~Vy + tj* — t>exp 


(30 


The evolution of the distribution functions f(x,v) and Fb(x, v) are then described by 
the Boltzmann equations for each distribution function. 


§1 

dt 


+ v^f + 


a • 


A / = §L X 

dv 1 St I coll 


(4a) 


m 

dt 


+ 


v •'$ Fb + a • 



6Fb 

St I coll 


(46) 
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In this simple model it is assumed that there are not forces acting on particles 
except during collisions (i.e. no plasma waves or plasma turbulence) — > a = 0, 
and we seek a stationary solution — ♦ d/dt = 0 

* • * f=f,\cotI <*«> 


V 


• ^ F b = 


6Fb 

St coll 


(5b) 


Physical interpretation: 

Particles travel along linear phase space trajectories except when removed 
from their respective distributions by a collision. Each collision removes a particle 
from the plasma contactor distribution and a particle from the ionospheric 
distribution 


&, = 

St I | coll 


( 6 ) 


(Only binary collisions are considered, and since the species chosen to represent the 
background distribution is a neutral atom, long-range interactions between particles 


are ignored.) Both 

ot 


coll 


and 


6 F b 
£t 


coll * n g eneral d e P en< i on x directly 


and indirectly through the value of the other distribution function. However, the 
sense in which we should consider the interaction of the distributions is weak is 


that we shall neglect variations in Fb in calculating CQ ^ . We can then 


use 


the value of / to calculate a solution for Fb and thus proceed iteratively until we 
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achieve some solution for / and F&. This iterative process may converge if the 
interaction of the two populations is sufficiently weak. The full iterative process 
will doubtless have to be carried out on a computer, but useful analytic results 
should be obtained by considering the first few iterations. The collisions of 
interest for this system are: 

1) atom-atom collisions 

2) atom-election collisions 

3) atom-ion collisions 

Since at least one of the particles is neutral in each of these three possibilities 
and the range of velocities is not very great, we will treat a collision (of whatever 
type) as being characterized by a collision cross section tr, independent of velocity 
(approximately) and different for each type of collisions. For the 1st iteration 
consider the density distribution in the plasma contactor cloud arising from Fk = 
const. The near free path of the plasma contactor cloud particles is 

A = 1 /No (7) 

where N is the number density of ionospheric particles. 

Estimate A numerically. For standard Shuttle conditions at an altitude of 
220 km ( N = 4xl0 10 cm' 3 ) these cross sections should be approximated by a 
typical gas kinetic cross section: o » 3xl0' 14 cm 2 , which gives 


A = (4xl0 1 °x3xl0- 14 )- 1 


cm 


833 cm 


8.33 m 
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which implies that we’re going to get a reasonable scale length for the cloud. 
For the conditions for TSS-1, the height is 300 km — ► N = 5xl0 8 cm' 3 — ► 


A = 667 m = 0.667 km 


the plasma contactor density in this limit may be written down by inspection: 


n(r) = no (r c /r) 2 e~ f / x 


( 8 ) 


We can calculate a characteristic scale size for the cloud from (8) 


n dr 



+ (n>/r ) 2 


(— l/A)e -f /* 


( 2 - M " 1 - - 1 


-r A 

rA 

V r \) 2/r + 1/A 


2A + r 

2A+ r 


note that as r— »oo, scale size— ► A; as r— >0, scale size— >r/2 


( 9 ) 


However, it is not surprising that the scale size is a function or r. 

We must now take this calculation to one higher iteration (at least) to get 
useful information about deviations from spherical symmetry 
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— > 

Fbo 



Initially consider that any scattering Which occurs between plasma contactor 
cloud particles and ionospheric particles has the effect of subtracting a particle 
from the ionospheric beam. Let 


Fbo = unperturbed ionospheric background distribution function 


F b 



= Fbo exp 


d x ' 


Ax>,y,z) i 


( 10 ) 


A(r) = 


1 

n(r)a 


with n(r)= initial iteration of 


contactor cloud distribution 


n(r) = n 0 ^r 0 /rj e r / A , where A = l/No 


is the initial mean free path estimated above. 


r = \/x 2 + y 2 + 2 2 , rt = \fx ' 2 + y 2 + z 2 
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Note that y,z = constant along an integration path. Let 6 = \f y 2 + z 2 , which 
has the role of an "impact parameter” 


Fb(^x,y,z^ = Fbo exp J dx > n(r')a^ 


= Fbo exp 


- norlo l* dx I 

J oo x ' 2 + b 2 


( 11 ) 


Begin by considering b = 0, i.e. the density distribution along the x axis. 


Fb(x) = Fbo exp 



exp(— x\/A) 


( 12 ) 


From standard integral tables, 



-1 e°* 

m-1 x m ~ l 


a f e 01 

m-1 J x” 1 - 1 


In this simple-minded model the density of the chemical release cloud is infinite at 
the origin, and so Fj(A) must go to zero at positive x. Will assume in what 
follows that when 6 = 0, x > 0. 



exp (~x[/A) 


(- 1 ) 


x< 



dx 


= - - + 

x 



i e~* dt 

c 
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= " l e ~ z/A + j 


(See Abramowitz and Stegun, Chapter 5 , pp. 227ff.) 
Ei(x) has a series expansion: 


(~) n x n 


as 


as 


Ei(x) = -7 - lnz - Y' i 

a ” ! 


x — ► 0 Ei(x)-+ oo 


x — *■ oo E\(x) — > 0 

Which guarantees proper limiting behavior for Fb(x) 


Fb{x) = Fbo exp J no 2 r 0 2 a 








(13) 


Note that (13) gives us Fb(x) < Fbo for all x < oo, (x > 0 ) as we would expect 
on physical grounds. 

Figure 1 displays the results of some calculations using the expressions 
derived above. The ionospheric particle density along the line of flight ahead of 
the plasma contactor (x measured in units of A) is shown for five different values 
of o = no 2 r 0 2 o/A: a= 10 , 3, 1 , 1/3, and 1/10. The curves for large a values lie 
below the curves for small s values as would be expected, since large a values 
correspond to larger values of the mass flux from the plasma contactor. 
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